TCHEBYCHEFF POLYNOMIALS AND THEIR 
RELATION TO CIRCULAR FUNCTIONS, 
BESSELFUNCTIONS AND LISSAJOUS-FIGURES 


by BALTH. VAN DER POL and TH. J. WEIJERS 


$1. Introduction 


In science it is usual to approximate an experimental function by 
a polynomial and then the question arises how to let the polynomi- 
al expression fit the observed data. To one of the many possible 
ways of solving this question attention was drawn by T cheb y- 
cheff1) anda Tchebycheff approximation is defined in 
such a way that, in the domain considered, the absolute value of the 
maximum deviation between the experimental curve and the ap- 
proximation (e.g. in the form of a polynomial of a given degree) 
should be as small as possible. Approximations in the T cheb y- 
cheff sense can obviously also be constructed to given analytical 
function. Tchebycheff himself investigated this approxi- 
mation in connection with the problem of the movement of the piston 
rod in an old Watt steam engine. The rod was connected to the beam 
by an arrangement of links in'such a way that the maximum devi- 
ation of the exact rectilinear motion should be a minimum ?). 

Obviously there are many modern technical problems where an 
application of an approximation in the Tchebycheff sense 
may be very useful. 

It is this approximation in ће Tchebycheff sense which is 
at the root of the definition ofthe Tchebycheff polynomials 
Т (и), because the latter is that polynomial of degree n, in which the 
coefficient of u^ is fixed and equal to unity, while theother coefficients 
are constructed in such a way that, in the domain —1 ux + l, 
the polynomial T,(u) is the best approximation in the T che b y- 


1) Tchebycheff, Œuvres, St. Petersburg 1899, page 111—143. 
2) 1. с. page 111—143, 533—538. 
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cheff sense to the function which is zero in the whole domain. 

Tchebycheff showed that it follows from the above con- 
ditions for the approximation to an arbitrary function that the 
absolute value of the maximum deviation between the given function 
and T, (и) occurs и + 1 times, the points p = + 1 included. 


In a general investigation into the best approximation to triode 
characteristics we were led to a study of Tchebycheff poly- 
nomials. It is the purpose of the present paper to expound their 
properties !) and to show the relation of these polynomials with 
circular functions, Bessel functions and Lissajous figures. 

In another paper we intend to describe the application of these 
polynomials to an approximate representation of triode character- 
istics. 


§ 2. The differential equation and its solution 


The Tchebycheff polynomials y = Т,(р) satisfy the differ- 
ential equation 


42у dy 
НО Speer iat di h= 
(1 equa tae ee 0. (1) 
Putting in this differential equation 
u = cos Ө, 
or 
0 = arc cos u, 
we obtain: 
d 
A nyo, © 
which is solved by 
yi etio \ 
2 = еі". j (3) 
Hence solutions of (1) are: 
yı = (ы Бім — 02)", 4 
i me Ee (4) 
ya = (и М1 — p?)”, | 


1) Inthe text books Tchebycheff polynomials are either not mentioned at all 
or they are treated in a very inadequate way. Short resumé's are to be found in C o u- 
rant-Hilbert, Mathematische Physik, Berlin (1926), and in Frank-Von Mises, 
Differentialgleichungen der Physik, Braunschweig (1930). 
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which could be taken as the two fundamental solutions of (1). In the 
same way, however, as it is sometimes more convenient to regard as 
fundamental solutions of (2) instead of (3) the expressions: 


= cos n 6, 


ya = sin n 0, 


we will in what follows consider as the two fundamental solutions 
of (1): 
ys = [+i V] — p" + (uy —i v1 — u3)"l 
1 (077—352 or 6 
э = z [p М 1—0)" — („— V1 — w^] | ©) 


or, which is the same thing: 


ys = cos (n arc cos y), 
ER | (7) 
ya = sin (n arc cos џи). 


Now it was shown by Tchebycheff that according to the 
definition of § 1 


уз = 277 T, (и). (8) 


The factor 2"—!, however, arises only from the fact that in the 
polynomial approximation to the function zero the coefficient of u” 
was chosen to be unity. As in, ће relation of Tchebycheff 
polynomials to other functions (such e.g. as already shown by (7)) it 
appears that this factor 2"—! can better be omitted, we will here use 
as the two basic solutions of (1): 


T, (p) =} [u + ?V/1— 2)" + (u—iv 1—y2)"] = cos (n arc cos y) 
and (9) 


О, (и) = 5 [(u.--24/1—p2)" — (иу 1—22)"] = sin (n avc cos y). 


In the present paper we will study the functions T,(u) and U,(u) 
only for integral values of м and in the real domain — 1 5 p x +1. 
The roots in (9) are defined as being positive, which obviates any 
ambiguity as to the multivaluedness of the function arc cos y. 

The definitions (9) at once show that 


T—, (и) = Ts (u), } 


U. (1) = — U. (ы). (19) 
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The definitions (9) suggest also a consideration of the functions 


V, (и) = cos (n avc sin |), \ (11) 


Ж, (и) = sin (n arc sin u), 
which are also solutions of (1), but it can easily be shown that 


V» (ш) = (— 1)” Ton (и), | 
Von+1 (ш) = (— 1)" Von +1 (ы), 
Wa (u) = (— )** Us (u), 

Wan: (и) = (— 1)” Тә +1 (и), 


so that we can limit ourselves to a consideration of the two functions 
Т, (и) and О, (и) as defined by (9) and from which it follows that 
T, (и) is a polynomial of degree n, while U, (u) is not. Both functions 
are real for integral values of n and for the domain of р considered, 
ie.—lzp,z-dl 

The definitions (9) also show that Т, and U, can be written as: 


n(n—3) 


(12) 


Т, (cos9) = 2"—! [cos" 9 — Hz 22 cos" 0 + 212 cos"—* 0 — 
n(n—4)(n—S) a 
— —3i28 — 005 $0-r...... ], 
ЕЕ (13) 
О, (cos 0) = 2"—! [cos"—! 9 — Te cos"—3 6 + 
A soss g MM) оет. ейн, 
or 
n n ut n(n—3) 5 
Т, (y) = iz t xax — 
n(n—4)(n—S) „_ 
— yt РТУ А 
n—2 (n—3) (n—4) v 
U, (и) = 2"—1 [u^ XA 1792 prs + EECIL ne pr-$— 
n—4 5) (n—6 — 
ВЕ. m Let... к=! 


Fig. 1 represents the first six T’s, fig. 2 the first six U’s, both as 
defined by (9), while figs. 3 and 4 give the same functions according 
to the old definitions, where each of the functions of figs. 1 and 2 are 
devided by 2"—1. 

Physica I 6 
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VIVA NAA ABE, 

TCH 
Wa UA 
UL NIAI 
CAN 


"A CR 
ДИА УАД 
LM DEUDA QAI 


-1 -0,8 -06 -04 -02 0, 
Fig.l. Tchebycheff polynomials with the normalisation used іп 
this paper. 


In the following tables we give the first ten T's and ten U's written 


out in full: 
Tou) — 1, Usu) —0, 
Tiy)— y. Up) = 1—2, 
Tau) 2293—1, Us(u) = 28у р, 
Talu) = 482—3џ, Uslu) = (402—1) /i— ps, 
T,(y) =8u'—8y?-+ 1, U (ш) = (885—441) /1—p2, 
Tilu) = 16—208 +54, Шш) = 169—129? 1) /1—p5, 
Т.) 23209 —48p5 + 182—1, Ош = (3205—32y5-- 64) /1—p, 
Т.) S 64u—112u54- 5653—7p., 0.40) —(640—80p*-F249—1) /1—pà— 
Talu) = 128—256 16011—32u? 4-1, Uau) —(12857—192454- 8095. —8u)4/1—tp5, 
Talu) —25699—576y + 432и5—120и5-Е9, U (u) = (2560—4480 + 24004—40р? + fie, — 


Tol) = 5121—1280 + 112005—40095 4-5042—1, (ш) = (512.°—1024 27 --67299— 160g 4- 102) V I— p" 
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| ү 
| PERY 


ІУ 
ИП ША ала аи 
“IC WA 
“WAC AAV 
MSNA 


Fig. 2. U-functions with the normalisation used in this paper 


It can further easily be shown with the aid of (1) that Т„ (u) and 
U, (и) can be written as 


n! 1—2 10% — n2\n h 
adn у же шыш 
Un (y) = (—1)—!.2°.я. 


nl ES 
(2n)! dui (1 — no^ 


Hence it follows that 
ат, (и) n 
= U, 
d Vic (u) | 
and (16) 
140, (ш) я 


du. V1 — 12 i 
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Fig.3. Tchebycheff polynomials with older normalisation. 


§ 3. Orthogonal properties 


The functions T, (u) and U, (и) further satisfy the orthogonal 
properties 


^41 du n 0, (mn) 
| Т, (u) Tn (0) a = cos й 0 cos n 0d 0 =. x | ) 
=! TE 0 25 mn 
a7 
+1 dy. T. | 0, (m zem) 
О, (u) Un (ы) —— = | sinn0sinm0d0 — т 
Jal V 1—p? 0 2 (m=n) 
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04 06 


Fig.4. U-functions with older normalisation. 


$4. Generating functions 
Also Т» (и) can be defined with the aid of a generating function: 


I — ty 


rape oe T. (u) PS (t< 1) (18) 


Another generating function, which is based on the older definition 
of Т» (и) (8) is 


1—72 


I—omcpR Bie =D) 


but this is based on the than exceptional value of To (и): 
To (p) = 1. 


Moreover, with our present definition to which we will further 
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adhere, it is easy to derive from (18) another generating function in 
the form of the logarithmic potential function as follows: 


eo i" 
lo _———— = У Т, E i«l 19 
Катора 6 (Ш... <D (19) 
By differentiating (19) with respect to u and using (16) we also 
find for the U's the generating function 
1 = 1 
= = b —— U, (p) t. (ф<1 20 
ini ew (<1) 00) 


$5. Recurrence relations 


From the definition (9) one can at once deduce the recurrence 
relations: 


Т»+а (и) — 2p Ts (и) + Ta (и) = 0; | (21) 


U, i (и) — 200, (и) + Usa (u) = 0. 


$6. Relation of the Tchebycheff polynomials T,(u) and the functions 
О, (и) with Besselfunctions 


We will now investigate the complete frequency spectrum (or the 
Fourier integral) of Besselfunctions. We therefore consider the 
general integral 

To 
Ku, m (о) = gon 2.0 dt, (22) 


—© 


where n — m is a positive integer. -The reduction of this integral is 
effected in the simplest manner with the aid of the symbolic calcu- 
lus ?). 

The method we use is the following. We temporarily introduce a 
variable x, which we shall eventually equate to unity, and which we 
represent symbolically as x = 1/5. We thereupon can easily eva- 
luate the integral as a function of ф. Further it is a simple matter 
then to write down the ,,original” as a function of x, which finally 
is replaced by x = }. 


1) Balth. van der Pol, Phil. Mag. 8, 861, 1929. 
Balth. van der Pol andK. Е. Міеѕѕеп, Phil. Mag. 8, 537, 1932. 
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We use the operational representation 


d 7, (2/5) = зз е 


1 
? 


or 
a ( xf m mea 
—— HJ, аыр» & 345, 
Е 2 2) \4ф 
where 
. 1 
Xx. 
Ф 


Thus the ,,image” of К is found as 


des. *) " а = 5 di 
n,m (о) = go 2( 2 Ф„\2 (222) i" +т а Du gen 2 (5) ә pcm d 


—©ф 


Putting further ¢ = 25v, we obtain 


+o : 1 
Ky, m (0) == gv t 2iov) v"—"dy = 
(о) | 25 Q2) 
. dn—» +o . 1 
— yan — piv? +2iwv) _ 7 
p f e VR dv 


or, with v + tw = u: 


K (о) —={"—т de DEEA еро? 7° ри du = 
nym * d m (25) = 

qm Дт" epu 

at VT qom (p) 


Now, according to the rules of the symbolic calculus ?): 


1 А х 


Р —%) 


(x — 02)" 


- and therefore 


— dis—3) ‚ (o? < x) 
0 (w? > x) 
The temporarily introduced variable x was equal to unity and 
hence 
+= _ OF [^7 f lev. em (1—«o2)—h, (02 <1) 
nn (0) = [ete ao 2H das (23 
Ma [ о (o? 1) 


1) loc.cit. 
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Using (15) the following special cases of (23) can now be written 
down: 


Ta (о) - 
Vine € os 
0 ; (> 1) 


\һ—1 U, (о) 2 
Kanı (о) -[ Tema л. [2 (— 1) E (02<1) (25) 
Е P ; (œ> 1) 


Р —4y'A/2x Р, (о), (w?<1 
К, +, En (e) - |. pem Таз Ше ШЫ £ Же | | | l ) (26) 
0 (о2> 1) 


2.(—1)" 


--9o 
Kno (a) = f em J, (t) dt = 


(24), (25) and (26) show that in the complete frequency spectrum of 
all Besselfunctions of integral order J,(!) and in general of #7, (0), 
where n — m is a positive integer, all frequencies greater than the 
angular frequency о = 1 are absent 1). And further that in the 
spectrum of J,() the frequencies between о and о + dw occur 
with an amplitude 


l T,(e) 
туло 


Similarly in the spectrum of t—'J,(¢) the frequencies between о 
and о + do occur with an amplitude 


1 U. (e) 


T п 


It further follows that, for о > + 1, the amplitudes of the com- 
ponents of J,(/ as given by (24) become infinite due to the long 
nearly sinusoidal part exhibited by the Besselfunction J,(/). If, 
however, we consider (25) and (26), which do not show the inifinity 
at о = + 1, we observe that the factor £—! in the expression 
1—17, (0) and £—* in the expression of J, ,;'(é) obviates the above 
resonance. 

The Fourier duality at once enables us to write the Bessel- 


1) We are aware that this property of J, (/) has also been noticed and connected with 
the physical properties of low-pass filters by Dr. Ph. 1e Corbeiller in his course 
of theoretical electricity at the Ecole supérieure des P.T.T. (Paris) 1933. 
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Fig. 5a. Fourier spectrum of J;(!) and Jj4(!). 
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functions as the integral over all the frequencies w inside the range 
—1<0< + 1. Thus: 


J,(é) = = 3)” | g tiot Te do, (27) 

Tut) (A fH usus) 

uec ие йш (28) 
ен = A | eto P, (o) do. (29) 


The amplitudes of the frequency spectra of I;(0), I,2(t), I7()/t, 


` 


L,2(t)/t, In, (£)] /t and I154,(0)/4/t are depicted in fig. 5, a, b and c. 


Fig. 5b. Fourier spectrum of 7—1J;(/) and £^lJ,^(!). 
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Fig. 5c. Fourier spectrum of /^ 7, (/) and 77.5, (t). 
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$7. Relations of Tchebycheff polynomials T,(u) and the functions 
U,,(p) with the Neumann- and Schlafli-polynomials 


Both the functions 7,(u) and U,(u) have a relation with the 
Neumann- and Schlafli-polynomials, which occur in the 
theory of Besselfunctions. 

The Neumann polynomial O,(z) 1) is defined by 


OF (2) = t; Te + мІ + x?)* + (x — V1 + e] e dx. 

With the aid of (9) this relation can be written: 

0, (2) = (Cf e Т, (ix) dx. 
Thus the symbolic representation follows: 

20, (5) = (~ i)" Т, (i9). (30) 
Similarly the Schläfli polynomial S,(z) ?) is defined by: 
2) = |, © + мІ + 2)" — (х — у + X?) e dx, 
0 V1 + x? 


and hence 


“2 —o tS @х. 
Sal) = 2-9 [oL 
Thus the symbolic representation follows: 


: U, (tx) 
S,(5) = 2 (— 171 ————. 31 
PS.) 2 (— 1757 (31) 
Further the known relation 5) 
1 С] 
1—4 = € Е, Jy, (2) O, (2), 
where & = 1, &, = 2 for n >O, 
can be written: 
РУ е, J,(2) £0, (9). (32) 
p —z n=0 


Both members of this expression can be regarded as the oper- 
ational representation of functions of x (where x = 1/2). 
1) Watson, Besselfunctions p. 278. | 


2 Watson, Besselfunctions p. 288. 
3) Watson, l.c. р. 271. 
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From the rules of the symbolic calculus it follows from (30) and 
(32): 


ё = У e, (— 0)" Ј, (2) T, (ix). 
0 


Substituting x = — iy, we obtain: 
e = Xe, (— 92" Jn (2) Т, (y) (33) 
n-0 


or, with the substitution y — cos 0, (33) becomes the fundamental 
generating function of Besselfunctions: 


giseosÜ — 


є, (— 2)" J, (2) cos ид, (34) 


ї [18 


n=0 


so that (32) can be regarded as the symbolic representation of the 
generating function (34) for Besselfunctions. 


§ 8. Relation of T,(u) and U,(u) with Lissajous figures 
The form of the definition (9) 


T,,(u) = cos (n arc cos u) 


is best suited to demonstrate the relation of this function with 
Lissajous figures. For let, fig. 6, a closed cosine curve with м periods 
be drawn on the surface of a glass cylinder with unit radius and with 
the z-axis as its axis. This special cosine curve is then given by 
y = 1, z = cos n Ө. Now consider the orthogonal projection of this 
curve on a plane parallel to the plane 0 = 0 through the axis of the 
cylinder. The abscis x of the projected figure is then x = p= cos 0 
and the ordinate remains z. Hence the projected figure gives the 
relation between џ = cos 0 and z = cos n 6 and is therefore repre- 


sented by 
z = cos (n arc cos p) = T,(u). 


It is to be remarked that, for the curve z = cos n 0 on the cylinder, 
the projection from the parts 0 < 0 < x and x <6 < 2r coincide. 
This is not the case when the figure on the cylinder is given by 
z — sin n 0, for then the projections from the front part, which, 
according to our definition (9), generates the U,-functions, have the 
opposite sign of z of those from the rear part, where x < arc cos i: 2m. 
This is all simply seen from a real glass model of the cylinder. 

When we further introduce the time 7 as variable instead of 0, the 
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projection of the curve z = cos n 0 on the cylinder can be written in 
the parametric from 

z = cos nt, 

X — cos Í, 
while the projection of the curve z = s?» n 0 on the cylinder can be 
written in the parametric form 

z = sin nt, 

X = cos t, 


and hence represent Lissajous figures for an integral ratio of 
the two composing frequencies. 


D 
^ 
ү 
1 
Lu 
"x 
[o 


D 


Fig. 6. Relation of Tchebycheff polynomial with Lissajous 
figure. 


From the above it follows, and this is clearly demonstrated by 
figs. 1 and 2, that on the one hand the Т,(џ) functions have a maxi- 
mum at the two sides u = + 1 and the curve going through а 
complete cycle retraces itself twice. On the other hand the curves for 
the U,(u) functions do not retrace themselves in a complete cycle 
but show in the second half of the cycle opposite values of those in the 
first half cycle. 
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$9. Fourier series and the functions T,(u) and U,(p) 


Suppose a periodic function z = f(t) be given. Then one complete 
period of this function can be imagined to be wound on a glass 
cylinder, whose radius is unity, so that the ends fit. Now write the 
function as z = /(0). Next consider the projection of this function 
again on the above mentioned plane parallel to the plane 8 = 0 
through the axis Z of the cylinder. In general a projection will thus 
be obtained, forming a closed loop (fig. 7), like the loops obtained 
on a cathode ray oscillograph, 
when a periodic function is 
studied with a cosinusoidal 
time-base. 

The part ABC may then be 
considered to be the projection 
of the domain 0 « 0< x of the 
original curve, and represented 
by z= 9,(u), and the part CDA 
is then tbe projection of the 
domain x < 0 < 2r of the origi- 
nal curve, and represented by Fig. 7. 

z = Q2 lu). 

The original function z = /(0) may be considered to be developped 

ina Fourier series: . 


z = У (a, cos n0 + b, sin n0). (35) 
n=0 
The projection of the individual components of (35) on the plane 


0 = 0 is derived from (35) by substituting » = cos 0. Hence we 
obtain on the plane 0 = 0: 


Ф:(и) = È [a Т„(и) + 5, U,u)); 


e) = 3 (a, Tale) — 5, Uu 


The curve AEC may be constructed by halving the ordinates of 
the loop. This curve is thus represented by 


(и) + ФА) ысу рдү. (36) 


n=0 
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and the coefficients of its development in a series of T che b y- 
cheff polynomials are equal to the coefficients of the cosine terms 
of the original Fourier series. 

The remaining top part of the loop, i.e. the function of p, for which 
the ordinates are given by the difference of the ordinates of the curve 
ABC and the middle curve AEC equals 


$6099 _ x b, Uu, (37) 
nol 


and the coefficients of its development in U,'s are equal to the 
coefficients of the sine terms in the original Fourier series. 

In experimental work sometimes the projected curve, i.e. the loop 
ABCDA is given and we wish to know the Fourier development 
of the ,original" function z = /(0). To this end we can develop thé 
middle curve AEC in the series (34) and the difference of ABC and 
AEC in the series (35). The coefficients a, and 5,, thus found, are 
the coefficients of the corresponding terms in the Fourier series. 

If the two branches ABC and ADC coincide, we have ф,(џ) =Ф(и), 
and the coefficients 5, vanish. 

Finally attention may be drawn to the fact that when a given 
function is represented by an infinite series of T,, polynomials and 
U, functions and when only say the first 2m terms of this series are 
considered, the approximation thus obtained is mot an approxi- 
mation in the Tchebycheff sense such as described in $ 1. 
Approximationsinthe Tchebycheff sense toa given function 
v(u) by a polynomial of. degree n can, according to the above 
considerations, be reduced to an approximation in the T che b y- 
cheff sense to a periodic function f(t) = ф (cos t) by a finite 
number of terms a,, cos m t and b, sin m t where the highest m 
occurring equals м. And as the wanted approximation is related to 
the absolute values of the maximum deviation between the approxi- 
mation and the given function, this would require a knowledge of 
the maxima and minima of finite Fourier series, about which 
subject not many data seem to be known yet. 
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